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The microquasar GRS 1758—258 is one of only three persistent black hole binaries that spend 
most of their time in the hard spectral state. It therefore provides the rare opportunity for an 
extensive long term study of this state. INTEGRAL has been monitoring the source since 2003 
spring during two ~ 3 months long Galactic Center viewing epochs each year, currently amount- 
ing to 1 1 epochs, including 2008 spring. Quasi-simultaneous RXTE observations are available as 
well. We present an analysis of the epoch averaged broad band spectra. Results include the ob- 
servation of clearly anti-correlated variations of the 3-300 keV flux and spectral softness, and of 
the occurrence a faint soft state, i.e., an extreme example of hysteretic behavior. The hard source 
spectrum and long exposures allow us to extend the analysis of the average hard state spectrum to 
~ 800 keV, resulting in the possible detection of a non-thermal Comptonization component. We 
also discuss the peculiar hardness-intensity diagram of this source. 
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1. Introduction 

GRS 1758—258 was discovered in 1990 during observations of the Galactic Center region 
performed with the Granat satellite [7, 16]. It is located only 0.66° from the bright low mass X-ray 
binary (LMXB) GX 5+1. Its average spectral and temporal properties — an exponentially cutoff 
power law spectrum (F= 1 .4—1.9, £f 0 id > 100 keV) and strong short term variability on frequencies 
up to 10 Hz, respectively — arc rather similar to those of Cygnus X-l [6, 4]. Based on these X- 
ray properties and on the detection of a double-lobed radio counterpart [10], GRS 1758—258 is 
considered to be a microquasar. It is one of only a handful persistent black hole binaries, with 
Cyg X-l and IE 1740.7—2942 being the only other ones predominantly in the hard state, making 
GRS 1758—258 a prime candidate for studying this complex accretion state. In contrast to the 
high mass X-ray binary (HMXB) Cyg X-l, however, GRS 1758—258 most likely has a low mass 
companion [11] and might thus belong to a source class that is usually transient due to Roche lobe 
overflow accretion. 

2. Long term lightcurves 

GRS 1758—258 has been monitored by RXTE from 1996 to 2007 with about two observations 
per week and by INTEGRAL from 2003 to today during two Galactic Center viewing epochs per 
year, each ~ 3 months long (Fig. 1, Fig. 2). Already the early RXTE monitoring revealed properties 
different from those of the canonical black hole hard state: in GRS 1758—258 the observed spectral 
hardness is correlated with the flux derivative on time scales of days (in the sense that the source 
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Figure 1: RXTE -PCA (observation-averaged) and 
INTEGRAL -ISGRI (3d binning) monitoring light 
curves of GRS 1758-258. 
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Figure 2: Close-up on the INTEGRAL -ISGRI 
light curve, including INTEGRAL -PICsIT detec- 
tions and epoch definitions (E01-E1 1). 


becomes softer with declining flux), while Cyg X-l shows no such hysteretic behavior [14]. The 
RXTE observations of 2001 also led to the discovery of the phenomenon of the rare faint soft state 
[12, 15]. Another occurrence of a faint soft state in 2003 spring is apparent in the light curves 
(Fig. 1, Fig. 2). Our earlier INTEGRAL-RXTE spectral fits already showed that, like the 2001 
faint soft state, this faint soft state is less luminous than the hard state, in this case by up to 
20% [8]. While this is again different from persistent FIMXBs, where the soft state is associated 
with higher bolometric luminosities, it is well within the range of hysteretic behavior displayed by 
LMXB transients (see, e.g., GX 339—4 [18]). In the following we report on updated results from 
modeling the epoch averaged broad band spectra as well as on using the hardest epochs to study 
the spectrum above 300 keV. 

3. Broad band spectra 
3.1 Epoch averaged spectra 

We have modeled the PCA and ISGRI/SPI/PICsIT broad band spectra averaged over each 
epoch 1 , successfully using either a cutoff power law or thermal Comptonization, plus a disk black 
body component when required, as well as a weak neutral iron Ka line. See [8] for an earlier 
version of similar broad band fits for the first four epochs. The PCA spectra are consistent with 
a column density given by interstellar absorption in the direction of GRS 1758—258, i.e., Ah = 
1.5 x 10 22 cm -2 . Typical exposures are ~ 30 ks for the PCA and ~ 1 —2 Ms for the INTEGRAL 
instruments (due to its smaller FOV JEM-X sees much less exposure from the source than the other 
INTEGRAL instruments and its data are not included here). Currently the fits for epochs 10 and 
epoch 1 1 still contain less data than those typical values since only data from the Galactic Centre 
Key Programme are included, more exposure will be available after the proprietary period. Note 
that epoch 10 was the last epoch covered by RXTE monitoring. 

1 ISGRI and SPI products were obtained with OS A 7.0, for PICsIT see [5], for the PCA see [6]. 
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Figure 3: Counts spectra and best fit models (upper panels) and best fit residuals (lower panels) for the 
quasi-simultaneous RXTE and INTEGRAL data of epoch 06 (left) and epoch 10 (right). The model in 
xspec notation is constant xphabs (diskbb+gaussian+powerxhighecut) . 


Fig. 3 shows the counts spectrum, best fit model, and residuals for two epochs, epoch 06, 
which shows a comparatively hard spectrum, and epoch 10, which is the second softest epoch after 
the faint soft state of 2003 (epoch 01). During these two softer epochs (01, 10) the source is not 
detected by PICsIT (for epoch 1 1 we do not have enough data yet). For the other epochs the most 
significant PICsIT detections, i.e., measurements in the energy band around ~300keV, are included 
in the fits. 



2003 2004 2005 2006 2007 2008 


Figure 4: Results of the cutoff power law fits for 
all epochs. Top: Fluxes (absorbed) in three energy 
bands. Middle: Flux constants for the PCA (red, set 
to 1), ISGRI (light blue, set to 1 in 2008), SPI (dark 
blue), and PICsIT (purple). Bottom : x^ cd values. 
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Figure 5: Results of the cutoff power law fits for all 
epochs. Top: Photon indices T (blue and red, the red 
cases indicate softer spectra) and 3-300 keV fluxes 
[arbitrary units]. Middle: Cutoff energies £cut- Bot- 
tom: Folding energies FfoW- 
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Figure 6: Long term average hard state spectrum of GRS 1758—258, obtained by averaging quasi- 
simultaneous RXTE- PCA and INTEGRAL -ISGRI/PICsIT data from the 2003-2007 epochs. Counts spectra 
and best fit models (upper panels) and best fit residuals (lower panels ) are shown for a thermal Comptoniza- 
tion fit without (left) and with (right) an additional power law component, i.e., in xspec notation the latter 
model is given by constant xphabs (diskbb+gaussian+compTT+power ) . 


Fig. 4 and Fig. 5 show results from the cutoff power law fits for all epochs. From Fig. 4 we 
see that the fluxes in several energy bands roughly follow the variations seen in the count rate light 
curves, that the cross calibration of the instruments generally agrees to within 10% (with the excep- 
tion of the PICsIT points which cannot be expected to constrain the normalization), and that x^ ed 
is < 1.5 (including 0.5% and 2% systematic uncertainties for the PCA and ISGRI, respectively). 
The upper panel of Fig. 5 shows that, on time scales of months, F and the 3-300 keV flux are 
anti-correlated with a correlation coefficient of —0.9. The middle panel and the lower panel show 
typical black hole hard state values for the cutoff and folding energies for epoch 02 to epoch 10. 
No cutoff can be detected in epoch 01 and epoch 11, where the latter is probably caused by the 
absence of PCA data and the low INTEGRAL exposure available so far. 

Results from thermal Comptonization fits for the first four epochs can be found in [8], Prelim- 
inary results for the remaining epochs show similar fit statistics and Comptonization parameters, 
i.e., describing a hot plasma with optical depths r ~ 0.5 — 1.5 and temperatures kT ~ 40 — 100 keV 

forced <!-5- 


3.2 Average hard state spectrum 

In the next step we averaged all PCA, ISGRI, and PICsIT spectra from epochs during which 
the source was clearly in the hard state, i.e., from epoch 02 to epoch 09. Currently PCA data 
sets are included starting at epoch 03 since a different PCU selection has been used for earlier data. 
(All PCA spectra from epoch 03 onward include only PCU 2 data, i.e., data from the best calibrated 
PCU. We are in the process of re-extracting the early epochs with this selection as well.) In addition 


5 


Monitoring GRS 1758—258 


Katja Pottschmidt 



channel energy (keV) 


Figure 7: Unfolded data and best fit model applying a cutoff power law plus an additional power law to the 
average hard state spectrum. The additional power law dominates below 10 keV and above 300 keV. 


the averaged SPI spectrum is still under investigation and is not yet included here. Caveats and 
ongoing work aside, we obtain a broad band spectrum with the most sensitive extension to energies 
> 300 keV to date, thanks to a significant source detection with PICsIT up to ~ 800 keV (for a 
detailed report on PICsIT observations, including those of GRS 1758—258, see [5]). 

Applying the two approaches used for modeling the individual epoch averaged spectra, the 
phenomenological cutoff power law and thermal Comptonization, to this long term average of the 
hard state spectrum, in both cases clearly cannot explain the data above ~ 300 keV. For the thermal 
Comptonization model we obtain X^A = 2-15 (Fig. 6, left) 2 . Another power law component was 
added to the models in order to describe these residuals. Interestingly this produced qualitatively 
different results for the cutoff power law and the thermal Comptonization model. 

In the latter case the best fit is a stable solution that reproduces the basic plasma parameters 
from the individual epoch fits — x = 1.4 ±0.1 and kT = 4 1 ^ keV for x^a = 1-43 — , adding a 
r = 1-4+Q4 power law that dominates the spectrum above 300 keV (Fig. 6, right). The presence 
of reflection in the spectrum can be constrained to being smaller than fl/2n < 20%. A weak 
~ 0.6 keV black body component is included in the model. This is a classical “hard tail” fit, 
indicating the possible presence of a non-thermal electron plasma, e.g., in the jet. We caution, 
however, that this is a preliminary result, there are calibration effects that still have to be evaluated, 
e.g., the question of the evolution of the ISGRI response with time. In addition the situation is less 
clear for the cutoff power law case: 

Here adding a power law leads to a local x} a \ minimum of 1.64, however, for a power law 
strength marginally consistent with zero and an unusual I of — 0.2 4 , 1 ) ' . On the other hand, a fit 
can be obtained that is formally as good as the thermal Comptonization plus weak power law fit, 
where the additional power law is strong, with an only roughly constrained I of 2.0 1 ] In this 
alternative fit the additional power law dominates below lOkeV and above 300 keV. The situation 

2 Systematic uncertainties for the PCA were increased to 1% in the long term average fits. 
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Figure 8: Summary of observational properties of the X-ray and radio emission of black hole X-ray binary 
transients, and their evolution over an outburst. Here the outburst states are classified by their position in the 
hardness-intensity diagram, i.e., in the “q” (central panel). 


is illustrated in Fig. 7. No blackbody component is required in this case. Although the cutoff 
power law fit is purely empirical, this result points at two different possibilities for the composition 
of the hot plasma. In order to further investigate this question we will apply models for hybrid 
thermal/non-thermal plasmas, e.g., eqpair, in future analyses of the average hard state spectrum. 

4. A special link in the evolution through the states 

Outbursts of black hole transients are often characterized by then - evolution in the hardness- 
intensity diagram (HID), typically concentrating on the ~ 2 — 20 keV energy range. A few years 
ago it was realized that outbursts principally follow similar tracks, with a hard rising phase - 
dominated by emission from the corona and/or jet — and a soft declining phase, dominated by 
emission from the accretion disk — , i.e., displaying considerable hysteresis, and with transitional 
phases at the end of the rise and the decline, where the spectral hardness changes strongly while 
the X-ray intensity is comparatively stable. Overall the HID is thus roughly “q-shaped” [2, 9]. In 
addition the different phases of radio jet emission over the outburst have also been characterized in 
relation to the HID. This especially includes the existence of the “jet-line”, a region in the upper 
transitional HID branch, where discrete ejections of jet material are most likely, and on the soft 
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Figure 9: Hardness-intensity diagram for the 

RXTE - PCA monitoring since 2003 (individual 
RXTE observations are plotted). Some of the 
epochs are highlighted in color. The faint soft state 
of epoch 01 plus several pointings before and after 
are shown in red to indicate the peculiar p-shape of 
the (anti-clockwise) evolution. 



Figure 10: Hardness-intensity diagram for the IN- 
TEGRAL -ISGR1 monitoring since 2003 (Id time 
bins), i.e., for the hard energy range. Compared to 
the < 20keV PCA HID, the > 20keV ISGRI HID 
does not show the hysteresis of a soft and a hard 
branch. Softening, however, is still associated with 
lower intensities. 


side of which the jet is quenched (see [17] for a recent example of the significance of the jet-line 
in Cyg X-l). Fig. 8 summarizes our current best knowledge of observational properties associated 
with the evolution of an outburst through the different states and the HT D . 

Persistent black hole binaries can be expected to occupy a comparatively small area of the “q”. 
For Cyg X-l monitoring with RXTE has, e.g., confirmed that the source spends most of its time 
on the hard state branch above 1-2% Z^dd with occasional excursions onto the upper transitional 
branch [17], never clearly descending into the soft thermal branch (“failed state transition”, “bright 
soft state”). Fig. 9 shows the HID of GRS 1758—258 for the RXTE-PCA monitoring during 
the INTEGRAL epochs. While it does show a brightening/hard and a declining/soft branch as 
transients do, this HID is rather “p-shaped” due to the faint soft state of epoch 01. Following 
the evolution into and out-off the faint soft state (Fig. 9, red points), a pattern similar to a regular 
outburst — developing anti-clockwise, from the hard branch over the upper transitional branch, 
down the soft branch, and thro ugh the lower transitional branch back to the hard branch — is 
observed. Consistent with the conclusion from spectral fitting (see section 1 and [8]) this suggests 
that the faint soft states correspond to the thermally dominated phase observed during LMXB 
outburst decays. 

We also show the rarely discussed HID for energies above 20keV (Fig. 10). Since the thermal 
disk component does not reach into this energy band there is only one hardness branch. Neverthe- 
less the effect of softening in the fainter states is still apparent, indicating that the faint soft state is 
not entirely due to a change in the relative contribution of the soft and hard spectral components, 
but that the latter also changes in itself. 

Interestingly GRS 1758—258 most likely has a low mass companion, namely a K0 III giant 
in an 18.45 ±0.10 day binary orbit, and is thus thought to be accreting via Roche lobe overflow 
[13, 11]. The accretion process can therefore be expected to show similarities to LMXB transients 
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and especially to include phases during which the mass transfer is reduced or has stopped com- 
pletely [14]. In any case GRS 1758—258 represents an important connection between persistent 
and transient sources 3 and the faint soft states provide a rare opportunity for studying what happens 
when the mass transfer is shut off, independent of the possible influence of a preceding outburst. 

5. Conclusions 

5.1 Summary 

The INTEGRAL and RXTE monitoring campaigns of GRS 1758—258 provide an exceptional 
database for studying the hard state of black holes. Above we presented the analysis of data from 
2003 to 2008, modeling spectra averaged over INTEGRAL viewing epochs as well as a long 
term average hard state spectrum, and discussing the hardness-intensity diagram above and below 
20 keV, thereby extending the study down to time scales of days. We report the following main 
results: 

• The epoch averaged 3-300 keV hard state spectra (typical INTEGRAL exposures of 1-2 Ms) 
can be modeled with a cutoff power law model or thermal Comptonization. They arc con- 
sistent with interstellar absorption. A weak iron line is detected and sometimes a weak soft 
excess. This picture is consistent with earlier results. 

• On time scales of months T, i.e., the softness of the broad band spectrum, and the 3-300 keV 
flux arc anti-correlated. 

• As reported in an earlier analysis of the first four epochs [8] , a faint soft state was observed 
in 2003 spring (epoch 01). Another softer-than-average state was observed in 2007 fall 
(epoch 10). 

• The 2003-2007 average hard state spectrum allows us to detect the source up to ~800keV 
with PICsIT. There are indications of a hard tail, i.e., the presence of a non-thermal Comp- 
tonization component, above 300 keV. In Cyg X- 1 such a component has been observed in 
softer states, i.e., not in the hard state (see, e.g., [1, 3]). 

• Although the HID of GRS 1758—258 is rather “p”- instead of “q”-shaped due to the faint soft 
state, the source can probably be understood in the canonical black hole transient picture, by 
interpreting the faint soft states as times of reduced mass accretion similar to outburst decays. 
This is consistent with the K-giant companion candidate, i.e., Roche-lobe overflow accretion. 

5.2 Outlook 

The analysis of this database is work in progress. Plans for future projects include to: 

• Refine the broad band spectra as more data and updated calibration becomes available. ( IN- 
TEGRAL monitoring of GRS 1758—258 is ongoing and additional data from recent epochs 
will become public.) 

3 IE 1740.7—2942 might be the only similar case. 
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• Investigate whether a model that includes non-thermal Comptonization can describe the long 
term average hard state spectrum (e.g., eqpair). Include SPI data in this analysis. 

• Place further constraints on the soft spectral component in the hard state by (i) evaluating the 
contribution of INTEGRAL -JEM-X data and (ii) using Swift -XRT data to constrain future 
epoch averaged spectra (a Swift AO-5 proposal has been submitted). Should INTEGRAL 
observe another faint soft state we plan to ask for a Swift ToO, in order to follow the evolution 
of the soft component. 

• Compare the > 20keV HID of GRS 1758—258 and other black hole binaries, as the HID of 
these sources has generally not been extensively studied at these energies. Model broad band 
spectra on time scales shorter than the INTEGRAL viewing epochs. 
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